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Abstract 

We give a physical interpretation of the entries of the reflection A'-matrices of 
Baxter's eight-vertex model in terms of an Ising interaction at an open bound- 
ary. Although the model still defies an exact solution we nevertheless obtain 
the exact surface free energy from a crossing-unitarity relation. The singular 
part of the surface energy is described by the critical exponents = 2 — ^ 
and «! = 1 — where /i controls the strength of the four-spin interaction. 
These values reduce to the known Ising exponents at the decoupling point 
/X = 7r/2 and confirm the scaling relations Us = oth + v and ai = a;, — 1. 
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Our understanding of phase transitions and critical phenomena has been greatly en- 
hanced by the study of exactly solved lattice models in statistical mechanics Chief 
among these models is Baxter's eight-vertex model, which exhibits continuously varying 
critical exponents 0. Such exact results provide valuable insights into the key theoretical 
developments of universality, renormalisation and scaling. The eight-vertex model is equiva- 
lent (see Figs. 1 and 2) to two Ising models coupled together by four-spin interactions . 
From the singular part of the bulk free energy scales as ~ \tY^^ as t — 0. Here t 
vanishes linearly with T — Tc, where Tc is the critical temperature. The variable /i measures 
the strength of the four-spin interaction M via 

exp(2M) = tan(/i/2). (1) 

When /i = 7r/m, where m is an even integer, the critical behaviour is modified to ~ 
|t|'^//^log|t|. This is the case in the Ising limit, where = 7r/2. The critical exponent 
describing the divergence of the bulk specific heat, ~ as t 0, is given by 

Of, = 2 — 7r//i, with a^, = (logarithmic) for the Ising model. 

A significant test of the scaling relations between critical exponents was given by Johnson, 
Krinsky and McCoy [0, who derived the correlation length exponent ^ = ^ for the eight- 
vertex model. Together with Baxter's result for this confirmed the validity of the bulk 
scaling law «{, = 2 — 2v. However, the situation is not so satisfactory for the surface 

critical behaviour as the eight- vertex model has not been solved for open boundary 

conditions as in Fig.2. Whereas integrability in the bulk is governed by solutions of the Yang- 
Baxter equation [5,0,0 , integrability in the presence of boundaries is governed by solutions of 



both the Yang-Baxter and reflection equations |10,11|. i^'-matrices satisfying the reflection 



equations have been found for the eight- vertex model |T2|-pr5[], but the diagonalisation of 
the transfer matrix remains a formidable problem. Here we nevertheless derive two surface 
critical exponents of the eight-vertex model, allowing a direct test of the proposed scaling 



relations between bulk and surface critical exponents 

The relation between the bulk Boltzmann weights a, 6, c, d of the eight-vertex model and 



the Ising couplings K, L, M is depicted in Fig.l. These weights are given by M 



a{u) = pQ6i^{X)9i{u)6i{\ — u), 

d{u) = poei{x)ei{u)e,{x-u), 



(2) 



where po is a normahzation factor. Here 9i{u) and O^^u) are the elhptic theta functions, 
e, (u) = 2 q'/^ sinh ^ ^ fl - 2g2" cosh ^ + q^A {1 - q' 

n=l 

04{u) = ft (l - 2g'"-^ cosh ^ + g^"-^) (1 - 



n=l 



of nome q = exp(— tt/'/J), where / and /' are the half-period magnitudes. In the principal 
regime, < m < A, with < A < /' and < g < 1. Here g — 1 at criticality. In terms of 
the vertex weights, the bulk Ising couplings are given by 

(3) 

2 

(4) 
(5) 



exp(4M) 



ac 




2 


bd ^ 






ad 


'9i{X- 


-u) 




- 


-u) 


ah 


'e^ixy 


2 


cd 







In the Ising limit, M = when A = -^I' , with the spectral variable u controlling the 
anisotropy of the Ising couplings K and L. 

For the lattice orientation of Fig. 2, the integrable boundary vertex weights can be 
written down from the entries of the i^-matrix. Now for the eight- vertex model this reflection 
matrix is of the general form K^{u) = K{u; r„), with 

■Kn{u) Ki2{u) 

where [jl5| 



(6) 
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Kii = p,--- r, (7) 

K o2,, M'^u) {r [dl{u) + dl{u)] - 9f{u) + 9l{u)} 



Here is a normalization factor and ^, r], r are arbitrary parameters. In principle, these 
three parameters are related to the surface couplings. We argue that the variable ^ controls 
the strength of the Ising surface coupling Ks- Similar to the bulk case, we see from Fig. 1 
that exp{AKs) is given by a ratio of the boundary weights Vij, which in turn are related to 
the K-matrix elements ||17|| , 



exp(4K.j - ^^^^^^ - ^^^(^/2)ir2i(V2)- ^ ^ 



The particular choice r = and ^ = |/' simplifies to 



exp(4ir,) = 



1 2 



(12) 



Comparison with the bulk parametrisation @ implies that the further choice of = —1 
leads to = Kg, i.e. equal bulk and surface couplings in the Ising spin formulation. These 
particular values of rj and r can be chosen for all ^, since the surface coupling Kg can be 
clearly set to be independent of rj and r, which can be seen from the product riir22. 

The surface free energy can be obtained by applying the inversion relation method, which 
is known to give the correct bulk free energy of the eight-vertex model (see, e.g. §13.7 of 
Ref. [|ll). By using the fusion procedure, the transfer matrix of the eight- vertex model with 
boundaries described by i^^-matrices has recently been found to satisfy a group of functional 



relations [18,19|. Ignoring the finite-size corrections, which are of no relevance here, the 



relations give the desired crossing-unit arity relation for the transfer matrix eigenvalues ^0 

A{u)A{u + A) = uj+{u)iu4u)p^^{u). (13) 

The factor 



Oi(A-M)Oi(A + i.) 



is a bulk contribution whereas the product ujj^{u)uj_{u) is a surface contribution, with |I8|JT9 

uj+{u) = Kti{u)b{-2u + X)K:^2iu + A) + Kt2{u)d{-2u + \)Kt2{u + A) 

-K+{u)a{-2u + A)ir+(w + A) - ir2^2(w)c(-2M + \)K+{u + A), (15) 
uj^{u) = K^iiu + A)rf(2n + \)K2i{u) + fs:22(M + X)h{2u + A)irfi(n) 

+ A)c(2n + A)irn(M) - fsTial^ + A)a(2M + A)fs:2i(^)- (16) 

Here K^{u) is the transpose of K^{—u + A) with replaced by etc. 

The bulk and surface free energies must both satisfy the crossing-unitarity relation (|13[). 
The surface energy can be separated from the bulk energy. As we are only predominately 
interested here in the surface critical behaviour rather than the precise form of the surface 
energy, we consider only the diagonal elements of the K-matrix. These terms are sufficient 
to extract the critical exponents and physically we do not anticipate any change in the 
critical behaviour arising from the off-diagonal terms ||2l|. Define A;, = and A^ = /€<j, 



then the bulk and surface free energies per site are defined by fb{u) = —logKb{u) and 
fs{u) = — logKs(M). From (13)-(16) we have 

K},{u)Kb{u + A) = p{u) (17) 

for the bulk and 

gi(e- - u)ei{i, + u)e^{i+ - u)9^{^+ + u) e,{2\ - 2u)e,{2\ + 2u) 

..(«).> + A) = (18) 

for the surface. 



We obtain the solution of (|T8|) for Ks(m) by applying the inversion relation method [|l| 



Let us first recall the derivation of Kb{u) from (|T^. It is convenient to introduce the variables 

a: = exp (— 7rA/2/) and w = exp (— vr-u/J) . (19) 

To obtain fb{w) the argument is to assume that Kb{w) is analytic and nonzero in the annulus 
x"^ < w < 1, allowing the Laurent expansion of fb{w), 



log Kk{w)= £ CnW^. (20) 



n=— oo 



Inserting this expansion into the logarithm of both sides of (|17D and equating the coefficients 
of powers of w then gives 

^^^"^ = - S n(l + a;-)(l-g-) • 

This is the desired result, from which the critical behaviour in the limit g ^ 1 is extracted 
by use of the Poisson summation formula In terms of the variable fi = nX/I', where 
/' ^ 7r/2 as g ^ 1, it follows that fh ~ p'^^^ as p 0, with ft, ~ p'^^^ logp if 7r//i is an even 
integer. Here the conjugate nome p = exp(— 27r///') vanishes linearly with the deviation 
from criticality variable t [Q]. 

We obtain the surface free energy by solving (|l^) under the same analyticity assumptions 
as for the bulk case, together with the further assumption that Hsiw) is analytic and nonzero 
in the annulus x < y± < 1, where we have defined y± = exp {—'n-^±/2I) . In this way we 
arrive at the result 

~^ n(l + x''^)[l — g^"^) 

n(l + x4")(l-g2n) • ^ ) 

Applying the Poisson summation formula leads to a series for fs in powers of the nome p. 

The phenomenology of critical behaviour at a surface is well developed P|7|JTB|. In this 
case two surface critical exponents can be obtained from the surface free energy; one from 
the surface specific heat, Cs ~ and the other from the "local" specific heat in the 

boundary layer, Ci ~ 1^1 Here the corresponding surface internal energy is given by 

es [P) ^ + ei (p) , (23) 

where ei{p) is the first layer internal energy, 

ei(p) TTT • (24) 
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The related specific heats follow as 



Cs-^ and (25) 
op op 



These definitions follow from P,[7|,p^ with the identifications p ~ t and ^± ^ Kg. From (p2D 
we find that as p — > 

e,(p)~p^~^ and ei(p)~p^/^. (26) 

As for the bulk logarithmic factor appears if vr/// is an even integer, with 

es(p) ~ P^"^ logp and ei(p) ~ p'^/^logp. (27) 

This behaviour is to be compared at /x = 7r/2 with the known Ising results where the 
logarithmic factor is observed [p^-p5|. 

In summary, we have derived the exact critical surface exponents 

TT TT 

= 2 and ai = 1 (28) 

2/i jJL 

for the eight-vertex model. At /x = 7r/2, = 1 (log) and ai = —1 (log), in agreement with 
the Ising results [^^-|25|. Recalling the bulk exponents = 2 — ^ and = ^ [§] W6 are 
thus able to provide a significant confirmation of the scaling relations [l6i[7|,|l6 



Os = C(b + ^ and ai = af, — 1 (29) 

between bulk and surface critical exponents. The derivation of other surface exponents 
awaits the diagonalisation of the transfer matrix, which remains a formidable open problem. 
We have found that the surface free energy scales as fs ~ p'^/'^^^ as p — > 0. It is interesting 



to observe that this is in agreement with the scaling behaviour of the interfacial tension [|26 
However, these two quantities differ away from criticality. 

Finally we note that, in the same spirit as this work, the critical magnetic surface expo- 
nent Sg of the two-dimensional Ising model in a magnetic field MM should also be obtainable 



from the dilute A3 model |27| , |28[| , which is known to lie in the same universality class as the 



Ising model in a magnetic field However, unlike for the eight- vertex model, where we 
have been able to disentangle the critical behaviour from the known solution of the reflection 
equations, the reflection matrices for the dilute ^3 model are yet to be determined. 
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FIGURES 
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t ^ rii = S gJf.+Mi+Ml 
r2i = S e-^^-^i+^' 



FIG. 1. The bulk and surface vertex and Ising spin configurations and their corresponding 
Boltzmann weights. The nearest-neighbour bulk interactions K and L are in the vertical and 
horizontal directions, respectively. The four-spin interaction is denoted by M and the general 
nearest-neighbour surface interactions by Kg^M] and M| in the vertical, SW-NE and SE-NW 
directions, respectively. The constants A and B do not enter into the critical properties. 
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FIG. 2. The geometric relation between the eight-vertex model lattice (dotted lines) and the 
Ising model lattice (broken and solid lines) . The Ising lattice is divided into two sub- lattices (solid 
and open circles). 
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